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A  novel  methodology  for  the  diagnosis  of  an  operating  proton  exchange  membrane  fuel  cell  (PEMFC) 
is  proposed  in  this  manuscript.  The  cell  impedance  can  be  represented  as  the  sum  of  two  terms:  the 
membrane  resistance  and  a  second-order  transfer  function.  An  experimental  procedure  and  a  method  to 
estimate  the  parameters  of  this  impedance  model  from  the  experimental  data,  by  applying  parameter 
identification  techniques,  are  proposed.  The  experimental  procedure  consists  in  reading  the  dynamic 
response  of  the  cell  voltage  and  current  after  the  occurrence  of  small  changes  in  the  load  value.  This 
experimental  procedure,  intended  for  application  on  operating  cells,  does  not  significantly  interfere  with 
the  cell  operation.  It  has  been  designed  supposing  that  the  frequency  range  of  the  relevant  cell  phenomena 
is  1  Hz  to  5  kHz.  The  required  experimental  equipment  is  inexpensive  and  easily  portable.  Analytical 
relationships  have  been  derived  in  order  to  calculate  the  following  cell  electrochemical  parameters  from 
the  model  parameters:  the  diffusion  resistance,  the  charge  transfer  resistance,  the  diffusion-related  time 
constant,  the  membrane  resistance  and  the  double  layer  capacitance. 

©  2008  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Our  society  is  suffering  from  fossil  fuel  shortage.  Fossil  fuels  (i.e., 
coal,  oil,  and  natural  gas)  also  contribute  to  a  number  of  environ¬ 
mental  problems  during  their  extraction,  transportation,  and  use. 
As  an  alternative  for  the  devices  consuming  fossil  fuels,  fuel  cells  are 
one  of  the  most  promising  means  of  producing  energy  in  portable 
systems.  In  particular,  the  use  of  proton  exchange  membrane  fuel 
cells  (PEMFCs)  in  automotive  systems  is  a  very  active  and  promis¬ 
ing  research  field.  Fuel  cell  vehicles  offer  substantial  reductions  in 
greenhouse  gas  emissions  [  1  ]  and  higher  energy  efficiency  than 
internal  combustion  engine  vehicles  [2], 

Mathematical  models  have  been  proposed  in  order  to  under¬ 
stand  the  physical  and  chemical  phenomena  involved  in  the  PEMFC 
operation.  These  include  ID  [3-6],  2D  [7]  and  3D  [8,9]  static  mod¬ 
els;  and  ID  [  10,11  ],  and  3D  [  12]  dynamic  models.  The  water  flooding 
in  the  PEMFC  cathode  has  been  modeled  from  physical-chemical 
first  principles  [13],  These  detailed  models  contain  a  large  number 
of  adjustable  parameters  and  their  computation  time  is  large. 

PEMFC  models  with  low  computational  cost,  intended  for  use  in 
control  applications,  have  also  been  developed.  In  this  context,  Ran¬ 
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dies  electric  models  and  electrical  impedance  spectroscopy  (EIS) 
techniques  are  employed  to  study  the  fuel  cell  behavior  and  to  esti¬ 
mate  the  value  of  the  cell  model  parameters  [14-22],  However,  EIS 
equipment  is  too  expensive  [23]  and  bulky  for  use  on  continuous 
in-field  assessment  of  operating  portable  PEMFC  systems.  There¬ 
fore,  this  approach  to  PEMFC  diagnosis  is  not  suited  for  in-field 
control  applications,  which  require  of  easily  portable  and  inexpen¬ 
sive  diagnosis  equipment,  able  to  be  integrated  into  the  portable 
power  systems. 

Parameter  estimation  techniques  have  been  successfully 
applied  by  other  authors.  An  ARX  model  was  proposed  in  [24] 
for  describing  the  dynamic  response  of  solid  oxide  fuel  cells. 
Other  approach  to  PEMFC  modeling  was  proposed  in  [25]:  the 
dynamic  response  of  the  PEMFC  is  represented  by  a  trans¬ 
fer  function  matrix  whose  inputs  are  the  cell  current  and  the 
oxygen  stoichiometry.  The  use  of  ARMAX  models,  with  the 
oxygen  and  hydrogen  flows  as  model  inputs,  and  PRBS  sig¬ 
nals  was  proposed  in  [26],  However,  these  approaches  neither 
provides  a  procedure  for  estimating  the  cell  physical-chemical 
parameters  from  the  cell  model,  nor  relates  the  estimated  transfer 
function  and  the  cell  impedance. 

A  method  for  diagnosis  of  PEMFC  through  current  interrup¬ 
tion  (Cl),  intended  to  be  used  in  portable  control  systems,  was 
proposed  in  [27],  It  exhibits  the  following  advantages  [27]:  the 
equipment  required  to  make  the  measurements  is  easily  portable 
and  inexpensive,  the  proposed  cell  equivalent-circuit  model  has 
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Nomenclature 

Qi  double  layer  capacitance  (F) 

Cg  concentration  in  cathode  active  layer  (mol  m-3 ) 

D  diffusion  coefficient  (m2  s-1 ) 

Foe  open  circuit  voltage  (V) 

f  Faraday  constant  (F) 

N  cycle  length  (bit) 

R  perfect  gas  constant  (J  mol-1  K-1 ) 

Rd  diffusion  resistance  (£2) 

Rm  membrane  resistance  (£2) 

Rp  charge  transfer  resistance  (£2) 

s  Laplace  transform  variable 

S  active  area  (m2) 

T  temperature  (K) 

Tsw  inverse  of  the  signal  rate  (s  bit-1 ) 

Zw  Warburg  impedance 

Greek  letters 

a  signal  amplitude  (V) 

S  diffusion  layer  thickness  (m) 

rja  anode  over-voltage  (V) 

rjc  cathode  over-voltage  (V) 

rjm  membrane  over-voltage  (V) 

rd  diffusion-related  time  constant  (s) 

c o  frequency  (rad  s-1 ) 


low  computational  cost,  and  the  PEMFC  electrochemical  param¬ 
eters  can  be  calculated  from  the  model  parameters.  However, 
this  diagnosis  method  exhibits  the  following  limitation.  As  the 
data  acquisition  method  for  parameter  estimation  requires  to 
completely  interrupt  the  cell  current,  this  cell  diagnosis  method 
interferes  with  the  cell  operation. 

The  main  objective  of  the  work  presented  in  this  manuscript  is 
to  propose  a  cell  diagnosis  methodology  with  all  the  advantages 
of  the  method  proposed  in  [27]  and  without  the  disadvantage  of 
interfering  with  the  cell  operation. 

A  simplification  of  the  PEMFC  model  presented  in  [27]  is  pro¬ 
posed  in  this  manuscript.  In  addition,  an  experimental  procedure 
and  a  method  to  estimate  the  parameters  of  this  simplified  model 
from  the  experimental  data,  by  applying  parameter  identification 
techniques,  are  proposed.  The  data  acquisition  procedure,  which 
can  be  applied  during  the  cell  operation,  consists  in  reading  the 
change  in  the  cell  voltage  and  current  induced  by  small  changes 
in  the  load  resistance.  Analytical  relationships  have  been  derived 
in  order  to  calculate  the  cell  electrochemical  parameters  from  the 
simplified  model  parameters.  The  proposed  diagnosis  methodol¬ 
ogy  allows  precise  estimation  of  the  following  cell  electrochemical 
parameters:  the  diffusion  resistance,  the  charge  transfer  resistance, 
the  diffusion-related  time  constant  and  the  membrane  resistance. 
The  model  simplification  is  obtained  at  the  cost  of  increasing  the 
error  in  the  estimation  of  the  double  layer  capacitance,  whose  value 
is  underestimated.  A  numeric  fitting  function  of  the  error  made  in 
the  estimation  of  the  double  layer  capacitance  is  provided.  This 
expression  allows  correcting  the  value  of  the  double  layer  capaci¬ 
tance,  reducing  significantly  the  estimation  error. 

2.  Experimental  set-up  and  procedure 

The  modelling  and  parametric  identification  of  the  PEMFC 
impedance  will  be  addressed  in  Section  3.  In  particular,  the  cell 
impedance  is  represented  by  Eq.  (1),  and  a  procedure  is  proposed 
for  estimating  the  parameters  Rm,  b,  c,  e  and/from  the  cell  experi¬ 


mental  data. 


Zcell(s)  =  Rm  + 


bs  +  c 
es2  +fs + 1 


(1) 


The  acquisition  of  the  experimental  data  is  discussed  in  this  sec¬ 
tion.  Essentially,  it  consists  in  reading  the  dynamic  response  of  the 
cell  voltage  and  current  after  the  occurrence  of  small  changes  in 
the  load  value.  This  experimental  procedure,  intended  for  applica¬ 
tion  on  operating  cells,  does  not  significantly  interfere  with  the  cell 
operation. 


2.1.  Experimental  set-up 

The  characteristics  of  the  PEMFC  used  in  the  study  are  100  cm2 
electrode  surface;  GORE  5761  membrane  with  membrane  thickness 
19  p,m;  SGL  10  BB  material  used  as  GDL;  and  serpentine  flow- 
field  at  the  anode  and  the  cathode  side.  The  operating  conditions 
are  the  cell  temperature  is  not  controlled  during  the  experiment; 
the  PEMFC  is  fed  with  air  and  hydrogen  at  atmospheric  pressure; 
the  hydrogen  flow  is  0.05  ml  min-1 ;  and  three  air  flows  are  used: 
0.61  ml  min-1, 0.51  ml  min-1  and  0.32  ml  min-1.  The  EIS  measure¬ 
ments  are  made  using  a  IM6  (ZAHNER  electrik)  workstation. 

The  PEMFC,  the  load  and  the  circuit  intended  to  produce  small 
changes  in  the  load  value  are  schematically  represented  in  Fig.  1  a. 
The  PEMFC  equivalent  circuit  is  the  subcircuit  named  ‘FC’.  The  load 
resistance  is  RL.  It  sets  the  cell  operating  point.  A  resistor  Rj  (with 
R/  »  Rl)  and  a  switch  represent  the  circuit  added  to  perform  the 
experimental  tests.  It  allows  forcing  small  changes  in  the  load  value, 
which  produce  small  excursions  from  the  cell  operating  point. 

An  implementation  of  this  test  bench  is  shown  in  Fig.  lb.  The 
switch  is  implemented  using  a  N-MOSFET,  whose  gate-to-source 
voltage  is  set  by  a  signal  generator.  As  the  N-MOSFET  channel  resis¬ 
tance  can  be  controlled  by  manipulating  the  gate-to-source  voltage, 
in  some  cases  the  resistor  R\  does  not  need  to  be  included  in  the 
circuit. 

A  microcontroller  with  digital  output  is  used  to  generate 
the  N-MOSFET  gate-to-source  voltage  (see  Fig.  lb).  The  voltage 
waveforms  used  during  the  identification  process  are  rectangu¬ 
lar  pulse-trains  and  pseudo-random  binary  sequences  (PRBSs).  The 
PRBS  signals  are  generated  by  using  a  shift  register  [28],  The  rea¬ 
sons  behind  the  selection  of  this  voltage  signals  will  be  discussed 
in  Sections  3.4  and  3.5. 

The  analyzer  device  can  be  designed  to  perform  the  paramet¬ 
ric  identification  calculations  described  in  Section  3.6.  Recursive 
least  squares  estimation,  widely  used  in  adaptive  control  [26,29], 
can  be  used.  Nevertheless,  the  results  discussed  in  this  manuscript 
have  been  obtained  by  using  off-line  identification,  which  requires 
storing  the  signals  in  the  memory  of  the  analyzer  device. 


2.2.  Experimental  procedure 

The  following  two  considerations  have  been  taken  into  account 
when  designing  the  experimental  procedure.  The  value  of  the  drain 
current,  while  the  N-MOSFET  is  in  the  on-state,  needs  to  be  small 
enough  to  justify  the  linear  approximation,  while  maintaining  an 
acceptable  signal-to-noise  ratio.  The  changes  in  the  cell  voltage 
induced  by  the  load  changes  should  be  smaller  than  the  thermal 
voltage  (i.e.,  approximately  26  mV  at  300  K)  [30], 

The  data  acquisition  procedure  consists  of  the  following  three 
steps: 

(1 )  Firstly,  a  5-kHz  rectangular  pulse-train  with  a  50%  duty  cycle  is 
applied  to  the  N-MOSFET  gate  for  80  ms.  The  cell  voltage  and 
current  is  read  with  a  sampling  period  of  10  pus.  This  data  is 
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Fig.  1.  Experimental  set-up:  (a)  schematic  representation  and  (b)  implementation  using  an  N-MOSFET  device. 


used  to  estimate  the  value  of  the  membrane  resistance  (Rm). 
The  procedure  is  discussed  in  Section  3.4. 

(2)  Next,  the  high-frequency  PRBS  voltage  waveform  described  in 
Table  1  is  applied  to  the  N-MOSFET  gate  for  0.15  s.  The  cell 
voltage  and  current  is  read  with  a  sampling  period  of  10  p,s. 

(3)  Finally,  the  low-frequency  PRBS  voltage  waveform  described  in 
Table  1  is  applied  to  the  N-MOSFET  gate  for  13  s.  The  cell  voltage 
and  current  is  read  with  a  sampling  period  of  500  jxs. 

The  data  obtained  from  Steps  (2)  and  (3)  is  used  to  estimate  the 
parameters  b,  c,  e  and  f  by  applying  the  parametric  identification 
techniques  described  in  Section  3.6. 

3.  Analysis  and  modeling 

3.1.  Main  assumptions 


(1)  The  oxygen  reduction  reaction  in  the  cathode  is  very  slow 
in  comparison  with  the  hydrogen  oxidation  reaction  [34,35]. 
Therefore,  the  anode  over-voltage  is  very  small  in  comparison 
with  the  cathode  over-voltage.  As  a  consequence,  the  anode 
over-voltage  contribution  to  the  cell  voltage  can  be  neglected 
in  the  model. 

(2)  The  double  layer  capacity,  which  is  usually  represented  by 
constant-phase  elements,  is  represented  by  a  pure,  single¬ 
frequency  theoretical  capacity  [16,22,33,36],  This  approxima¬ 
tion,  which  reduces  the  computational  cost  of  the  model,  is 
reasonably  accurate  at  low  and  medium  frequencies  (i.e.,  fre¬ 
quencies  below  a  few  hundred  Hz).  The  error  due  to  this 
hypothesis  is  justified  for  the  sake  of  obtaining  a  model  simple 
enough  to  be  suited  for  control  applications. 

3.2.  Modelling  of  the  Warburg  impedance 


The  voltage  drop  across  the  fuel  cell  (Vceu)  can  be  written  as  a 
function  of  the  steady-state  open-circuit  voltage  of  the  cell  (E0c). 
and  the  over- voltages  of  the  anode  {rjf),  the  cathode  (pc)  and  the 
membrane  (rjm)  [31,32]: 

VCell  =  £oc  —  l?»-  lta-9c  (2) 

The  open-circuit  voltage  (Eoc )  is  modeled  as  an  ideal  voltage  source. 
The  external  operating  conditions  and  the  cell  current  are  consid¬ 
ered  constant  during  the  experiments.  Therefore,  the  over-voltages 
(p3,  Pm,  pc)  are  assumed  to  be  only  dependent  on  the  FC  inter¬ 
nal  phenomena.  The  over-voltages  of  the  different  layers  can  be 
modeled  as  electric  components  that  represent  the  FC  dynamic 
behavior. 

The  cell  voltage  in  Eq.  (2)  can  be  modeled  by  the  electric  circuit 
shown  in  Fig.  2a,  which  is  composed  of  Randles  models  connected 
in  series  [15,33  ].Z^  and  are  the  Warburg  impedances  associated 
to  the  gas  diffusion  in  the  anode  and  the  cathode,  respectively.  Rp 
and  Rp  are  the  charge  transfer  resistances  in  the  anode  and  the 
cathode.  and  are  the  double  layer  capacities  in  the  anode 
and  the  cathode.  Finally,  Rm  is  the  membrane  resistance. 

Two  additional  hypotheses  are  made  in  order  to  simplify  the  cell 
model  shown  in  Fig.  2a.  As  a  result,  the  simplified  model  shown  in 
Fig.  2b  is  obtained. 


The  Warburg  impedance  (Zw)  can  be  written  in  the  Laplace 
domain  as  a  function  of  the  finite  length  diffusion  [15]: 

Z„W-X,tan^S  (3) 

Vs7rd 

where  the  diffusion  resistance  (Rd )  and  the  diffusion  time  constant 
(rd)  can  be  calculated  from  the  following  expressions: 


RTS 

d  SCgDn2F2 


(4) 


(5) 


The  approximation  of  the  Warburg  impedance  shown  in  Eq.  (6)  was 
proposed  in  [27],  This  approximation  is  equivalent  to  model  the  cell 
by  using  the  circuit  shown  in  Fig.  3.  The  value  of  the  parameters  Ri , 
R2,  Ci  and  C2  were  calculated  in  [27]  and  they  are  shown  in  Table  2. 


Zw(s)  Rd  ( 1  +  RiCitdS  +  1  +  R2C2rds) 


(6) 


The  exact  value  of  the  impedance,  calculated  from  Eq.  (3),  and  the 
approximated  value,  obtained  from  Eq.  (6)  using  the  previously  cal¬ 
culated  values  of R\,R2,  Q  and  C2,  are  plotted  in  Fig.  4.  The  absolute 
error  is  also  shown  in  Fig.  4.  Frequencies  (<w)  in  the  range  from  1 
rad  s_1  to  104rads_1and  diffusion-related  time  constants  (rd)  in 


Table  1 

Parameters  defining  the  two  PRBS  voltage  waveforms  applied  during  the  data  acqui-  Table  2 

sition  process  Fitted  values  of  the  parameters  in  Eq.  (6) 


Units 


Unitless 

Unitless 

Unitless 

Unitless 
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Fig.  3.  Equivalent  circuit  model  of  the  PEMFC  with  the  Warburg  impedance  calcu¬ 
lated  from  Eq.  (6). 


the  range  from  0.2  s  to  1  s  have  been  considered.  As  a  consequence, 
the  range  of  jwrd  in  Fig.  4  is  0.2  to  104. 


3.3.  Modelling  of  the  PEMFC  impedance 


The  cell  impedance  (Zcell),  calculated  from  the  circuit  shown  in 
Fig.  3,  is  as  follows: 


Zcen(s)  =  Rm  +  - 


1 


sQu  +  (l/^p  +  ^w) 
It  can  be  written  as  follows: 

as2  +  bs  +  c 


Zcell(s)  =  Rm  + 


ds?+eS2+fi+g 


where  the  parameters  a,  b,  c,  d,  e  and /are  related  to  the  cell  elec¬ 
trochemical  parameters  as  shown  below  [27], 


a  =  9.76  x  10_4Rprd  (9) 

b  =  0.304Rprd  +  3.38  x  10-2Rdtd  (10) 

c  =  Rp  +  0.949Rd  (11) 

d  =  9.76  x  10-4CdiRptj  (12) 


Fig.  4.  Warburg  impedance  calculated  from  Eq.  (3)(-)  and  Eq.  (6)(  x ).  Absolute  error. 


e  =  3.37  x  10-2CdlRdTd  +  0.3048CdlRpTd  +  9.76  x  W~4rj  (13) 
/  =  CdlRp  +  0.3048rd  +  0.949CdlRd  (14) 

g  =  l  (15) 

A  further  approximation  for  the  cell  impedance  is  proposed  in 
this  manuscript,  which  reduces  the  order  of  the  model  and,  con¬ 
sequently,  its  computation  time.  The  parameters  a  and  d,  which  are 
two  orders  of  magnitude  smaller  than  the  other  parameters  (i.e.,  b, 
c,  e,f  and  g),  are  neglected  in  Eq.  (8).  Making  this  approximation, 
the  following  model  for  the  cell  impedance  is  obtained: 

Zcell(s)  =  Rm+ (16) 

Once  the  parameters  b,  c,  e  and/have  been  estimated  from  the  cell 
experimental  data,  the  following  cell  electrochemical  parameters 
can  be  calculated  from  Eqs.  (10),  (11 ),  (13)  and  (14):  the  double  layer 
capacitance  (Cdi),  the  diffusion  resistance  (Rd),  the  charge  transfer 
resistance  (Rp)  and  the  diffusion-related  time  constant  (rd). 

The  models  described  by  Eqs.  (8)  and  (16)  predict  the  same 
steady-state  value  of  the  cell  impedance: 

lim[Zcell(s)]  =Rm  +  c  =  Rm+R„  +  0.949 Rd  (17) 

S-r0 

This  value  of  the  steady-state  impedance  can  be  calculated  directly 
from  the  circuit  shown  in  Fig.  3,  by  replacing  the  capacitors  by  open- 
circuits.  As Ri  =  0.8463  andR2  =  0.1033  [27],  the  value  obtained  by 
adding  R,  and  R2  is  0.949. 

The  proposed  simplification  of  the  impedance  model  does 
not  introduce  additional  error  in  the  estimation  of  the  diffusion 
resistance  (Rd ),  the  charge  transfer  resistance  (Rp )  and  the  diffusion- 
related  time  constant  (rd).  The  model  simplification  is  made  at  the 
cost  of  increasing  the  error  in  the  estimation  of  the  double  layer 
capacitance  (Cdi),  whose  value  is  underestimated. 

The  error  in  the  estimation  of  the  double  layer  capacitance  (eCdl ) 
is  approximately  described  by  Eq.  (18).  Therefore,  the  double  layer 
capacitance  is  estimated  by  adding  the  following  two  terms :  (a)  the 
value  calculated  from  the  estimated  model  parameters  (i.e.,  b,  c,  e 
and/)  using  Eqs.  ( 10),  ( 11 ),  ( 13)  and  ( 14)  and  (b)  the  estimation  error 
calculated  from  the  following  equation: 

£cdl  =  -0.005155Rp 0  9197 rd  -  0.2629Rdrd  +  0.02855rd  (18) 

3.4.  Estimation  of  the  membrane  resistance 

As  it  was  described  in  Section  2.2,  the  experimental  data 
required  to  estimate  the  membrane  resistance  (Rm)  is  obtained  by 
applying  a  5  kHz  rectangular  pulse-train  voltage  to  the  N-MOSFET 
gate  for  80  ms.  The  cell  current  and  voltage  are  read  with  a  sampling 
period  of  10  pus. 

The  Rm  term  in  Eq.  (16)  constitutes  a  proportionality  factor 
between  the  changes  in  the  cell  current  and  voltage.  The  effect 
of  the  membrane  resistance  manifests  at  high  frequencies,  usually 
overlapping  with  the  high-frequency  inductive  behavior,  which  is 
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13.5 


AI 


time  (ms) 


Fig.  5.  Simulated  change  in  the  cell  current  (AI)  and  voltage  (AV)  in  response  to 
the  applied  N-MOSFET  gate-to-source  voltage  (i.e.,  5  kHz  pulse-train  with  50%  duty 
cycle).  The  membrane  resistance  is  Rm  ~  Jpjp 


represented  by  Z,  in  Eq.  (19).  The  method  used  to  estimate  the  mem¬ 
brane  resistance  is  described  by  Eq.  (20).  An  example  is  shown  in 
Fig.  5. 


UmZceU(s)  =  Rm+Zi(s)  (19) 

r  i  |av| 

Re  |jimZcell(s)J  =  Rm  ~  ^1  (20) 


3.5.  Selection  of  the  PRBS  voltage  signals 


during  10  cycles  of  the  PRBS  voltage  waveform,  the  data  acquisition 
process  takes  at  least  10 TSWN.  For  this  frequency  range  selection 
(i.e.,®inf  =  2n  rads-1  -*•  TSWN  =  2jr/®jnf  =  1  s  cycle-1 ),  the  process 
should  take  at  least  10  s.  Considering  that  the  sampling  frequency 
needs  to  be  at  least  4.5  x  104  samples  s-1 ,  the  total  number  of  sam¬ 
ples  should  be  at  least  4.5  x  105  samples. 

An  alternative  approach  consists  in  first  applying  10  cycles  of  a 
low-frequency  PRBS  signal,  whose  frequency  range  is  0.77-75  Hz, 
and  next  applying  10  cycles  of  a  high-frequency  PRBS  signal,  whose 
frequency  range  is  71  Hz  to  2  kHz.  The  characteristics  of  these  two 
waveforms  are  summarized  in  Table  1  and  their  frequency  spec¬ 
trum  is  shown  in  Fig.  6.  This  approach  allows  reducing  the  number 
of  samples  in  an  order  of  magnitude  (i.e.,  the  required  number 
of  samples  is  now  approximately  4  x  104  samples).  As  a  conse¬ 
quence,  using  only  one  PRBS  signal  to  cover  all  the  frequency  range 
is  impractical,  because  it  leads  to  an  unnecessarily  large  number  of 
samples. 


3.6.  Estimation  of  the  parameters  b,  c,  e  andtf 


The  cell  voltage  and  current  measured  during  the  application  of 
the  two  PRBS  voltage  signals  are  used  to  estimate  the  impedance 
parameters  b,  c,  e  and/.  As  Zcen(s)  =  14e[](s)//cei[(s),  the  parametric 
identification  is  performed  by  considering  that  the  cell  current  (lcen ) 
is  the  process  input  and  the  cell  voltage  (Vcell)  is  the  process  output. 

The  effect  of  the  membrane  resistance  is  substracted  from  the 
measured  cell  voltage  (l/ceii)  as  shown  in  Eq.  (23).  The  cell  current 
(fceii)  and  the  calculated  cell  voltage  (V*ell)  are  used  to  estimate  the 
parameters  b,  c,  e  and  /  on  the  basis  of  Eq.  (24).  To  this  end,  the 
frequency  domain  system  identification  toolbox  for  Matlab,  called 
Ident,  has  been  used. 


Rm/Cell 


(23) 


Vcel,W  = 


;  ftg  +  c 
es’  -fs-ri 


fcell(s) 


(24) 


In  order  to  estimate  the  parameters  b,  c,  e  and/in  Eq.  (16),  two 
PRBS  voltage  waveforms  are  applied  to  the  N-MOSFET  gate  (see 
Section  2.2).  The  reason  for  using  two  PRBS  signals  instead  of  only 
one  PRBS  signal  is  reducing  the  number  of  collected  experimental 
samples.  In  order  to  justify  this  assertion,  the  use  of  only  one  PRBS 
voltage  waveform  is  considered  first. 

The  frequency  range  for  the  relevant  cell  phenomena  is  sup¬ 
posed  to  be  1  Hz  to  2  kHz.  On  the  other  hand,  the  spectral  power 
of  a  PRBS  signal  generated  using  an  shift  register  is  described  by 
Eq.  (21  )[28],  The  signal  frequency  range,  [a>inf,  &)sup],  is  shown  in 
Eq.  (22),  where  N  is  the  cycle  length,  Tsw  is  the  inverse  of  the  signal 
rate  and  a  is  the  signal  amplitude. 


^M-a2(N  +  1)Tswpin(ftirsw)12 

1  N  [  tuTSw/2  J 


2n 

TSWN 


2.8 

S  a)  <  — - <wSup 


(21) 

(22) 


Once  the  desired  frequency  range  has  been  selected  (i.e.,  the  values 
of  ft)jnf  and  <ySup  have  been  set),  the  corresponding  values  of  Tsw  and 
N  can  be  calculated  from  Eq.  (22).  In  order  to  cover  the  frequency 
range  of  the  cell  relevant  phenomena  (i.e.,  1  Hz  to  2  kHz),  the  fre¬ 
quency  range  of  the  PRBS  waveform  should  be  [winf  =  2n,  u>sup  = 
4;rl03]rad  s-1. 

While  the  PRBS  voltage  waveforms  are  applied  to  the  N-MOSFET 
gate,  the  cell  voltage  and  current  need  to  be  read  with  a  sampling 
frequency  of  at  least  10  times  Tsw.  As  a  consequence,  the  sampling 
frequency  needs  to  be  at  least  4.5  x  104sampless_1.  On  the  other 
hand,  adopting  the  heuristic  criterium  of  collecting  data  at  least 


The  Ident  toolbox  for  Matlab  allows  estimating  the  coefficients  of 
the  transfer  function  shown  in  the  following  equation: 


(25) 


The  parameters  b,  c,  e  and /can  be  calculated  from  the  coefficients 
k,  Tz,  Tpi,  Tp2  using  the  following  equations: 


c  =  k  (26) 

b  =  kTz  (27) 

e  =  TplTp2  (28) 

/  =  Tpi  +  Tp2  (29) 
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Fig.  9.  EIS  spectrum  corresponding  to  the  following  air  flows:  (A)0.61  mlmin  1 ,  (B) 
0.51  mlmin-1  and  (C)  0.32  ml  min-1 ;  and  fitted  third-order  models  (-  -). 


Table  3 

Cell  parameters  corresponding  to  successive  stages  in  the  cathode  flooding  process 


4.  Results  and  discussion 

The  test  circuit  shown  in  Fig.  lb  allows  obtaining  switching 
times  under  10-5  s.  An  example  is  shown  in  Fig.  7.  As  the  time 
constant  of  the  cell  response  is  typically  above  10-4  s,  switching 
can  be  considered  instantaneous. 

The  PRBS  voltage  applied  to  the  N-MOSFET  gate  induces  changes 
in  the  cell  current  whose  waveform  is  not  exactly  a  PRBS  (see  Fig.  8 
a).  Nevertheless,  the  cell  current  exhibits  a  frequency  response 
similar  to  the  frequency  response  of  a  PRBS  (see  Fig.  8b).  PRBS 
changes  in  the  cell  current  could  be  obtained  by  using  galvanostatic 
equipment.  However,  the  proposed  procedure  allows  obtaining  the 
desired  result  by  using  inexpensive  and  easy-to-use  equipment. 

In  order  to  compare  the  results  obtained  from  the  third-order 
and  the  second-order  impedance  models,  described  by  Eqs.  (8)  and 
(16),  respectively,  the  following  study  has  been  performed. 

EIS  measurements  have  been  made  on  the  cell  and  under  the 
experimental  conditions  described  in  Section  2.1 .  In  order  to  induce 
different  stages  of  the  flooding  process,  three  air  flows  have  been 
used.  The  results  are  shown  in  Fig.  9.  The  third-order  impedance 
model  has  been  fitted  to  the  experimental  data  corresponding  to 
these  flooding  process  stages,  named  A,  B  and  C.  The  electrochem¬ 
ical  parameters  of  the  cell,  calculated  from  the  fitted  models,  are 
shown  in  Table  3. 


1  jl  y  5  KHz  (square  signal) 


-0  5  PRBS  (high)  PRBS  (low) 

2  6  10 
time(s) 


Fig.  10.  State  of  the  switch  that  controls  the  load  changes  during  the  data  acquisition 
process  (1:  closed;  0:  open). 


Firstly,  the  data  acquisition  procedure  and  the  cell  voltage  have 
been  simulated  at  these  three  different  stages  of  the  cathode  flood¬ 
ing  process,  using  the  third-order  impedance  model.  The  state  of 
the  switch  that  controls  the  load  changes  during  the  simulated 
experimental  run  is  represented  in  Fig.  10.  The  pulse-train  voltage 
waveform  is  applied  to  the  N-MOSFET  gate  in  the  first  place.  Next, 
the  two  PRBS  voltage  waveforms  are  applied. 

Secondly,  uniformly  distributed  noise  has  been  added  to  the  cell 
voltage  and  current  values  obtained  by  simulating  the  experimen¬ 
tal  procedure.  Each  simulated  value  of  the  cell  current  has  been 
incremented  by  a  U[- 50,  50)  mA  distributed  random  variate.  Anal¬ 
ogously,  each  simulated  cell  voltage  has  been  incremented  by  a 
U(— 1,  l)mV  distributed  random  variate. 

For  instance,  the  cell  voltage  calculated  from  the  third-order 
impedance  model  during  the  application  of  the  low-frequency  PRBS 
voltage,  and  modified  by  adding  the  corresponding  l/(-l,  l)mV 
distributed  random  variate,  is  represented  by  dots  in  Fig.  11.  The 
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5.  Conclusions 


Fig.  11.  Cell  voltage  during  the  application  of  the  low-frequency  PRBS  voltage  wave¬ 
form.  Simulated  experimental  data  (dots)  corresponding  to  the  “A"  flooding  process 
stage  and  values  predicted  by  using  the  fitted  second-order  impedance  model  (con- 


Table  4 

Parameters  estimated  using  the  diagnosis  method  proposed  in  this  manuscript 


impedance  model  parameters  correspond  to  the  A  stage  of  the 
cathode  flooding  process  (see  Table  3). 

Thirdly,  the  cell  voltage  and  current  values  previously  calculated 
have  been  used  to  estimate  the  parameters  of  the  second-order 
impedance  model.  Three  impedance  models  are  obtained,  corre¬ 
sponding  to  the  three  flooding  process  stages.  The  Nyquist  and 
Bode  diagrams  for  the  third-order  and  the  second-order  impedance 
models  are  shown  in  Fig.  12.  Continuing  with  the  example  shown 
in  Fig.  11,  the  corresponding  cell  voltage  predicted  by  the  second- 
order  impedance  model  is  represented  by  a  continuous  line  curve 
in  Fig.  11. 

Finally,  the  cell  electrochemical  parameters  predicted  by  the 
second-order  impedance  models  are  calculated.  They  are  shown 
in  Table  4.  The  values  of  the  double  layer  capacitance  have  been 
corrected  by  adding  the  corresponding  estimated  errors,  which  are 
evaluated  from  Eq.  (18). 

Comparing  the  results  shown  in  Tables  3  and  4,  the  conclusion 
is  that  the  model  simplification  and  the  data  acquisition  procedure 
proposed  in  this  manuscript  do  not  introduce  additional  error  in 
the  estimation  of  the  following  electrochemical  parameters:  the 
membrane  resistance  (Rm),  the  diffusion  resistance  (Rj),  the  charge 
transfer  resistance  ( Rp )  and  the  diffusion-related  time  constant  (rq). 
The  impedance  model  simplification  is  made  at  the  cost  of  increas¬ 
ing  the  error  in  the  estimation  of  the  double  layer  capacitance  (Cji), 
whose  value  is  underestimated.  A  more  precise  estimation  of  the 
double  layer  capacitance  is  obtained  by  adding  to  the  estimated 
value  the  correction  term  described  by  Eq.  (18). 


Re(Z)  (Si.)  frequency  (red  s'1) 


A  method  for  electrochemical  parameter  estimation  in  operat¬ 
ing  PEMFC  has  been  discussed  in  this  manuscript.  It  has  the  three 
main  advantages  of  the  method  proposed  in  [27],  Firstly,  the  equip¬ 
ment  required  to  make  the  experimental  measurements  is  easily 
portable  and  inexpensive.  Secondly,  the  proposed  cell  impedance 
model  is  well-suited  for  control  applications:  it  has  low  compu¬ 
tational  cost  and  it  allows  reproducing  the  dynamic  behavior  of 
the  PEMFC.  Finally,  the  following  PEMFC  electrochemical  param¬ 
eters  can  be  calculated  from  the  model  parameters:  the  diffusion 
resistance,  the  charge  transfer  resistance,  the  diffusion-related  time 
constant,  the  membrane  resistance  and  the  double  layer  capaci¬ 
tance.  The  proposed  diagnosis  method  has  an  additional  advantage: 
the  experimental  measurements  can  be  performed  during  the  cell 
operation.  The  data  acquisition  procedure  does  not  interfere  with 
the  cell  operation.  It  consists  in  reading  the  change  in  the  cell 
voltage  and  current  induced  by  small  changes  in  the  load  resis¬ 
tance. 

Acknowledgments 

This  work  has  been  supported  by  the  Spanish  CICYT,  under  DPI- 
2007-61068  grant,  and  by  the  IV  PRICIT  (Plan  Regional  de  Ciencia 
y  Tecnologia  de  la  Comunidad  de  Madrid,  2005-2008),  under  S- 
0505/DPI/0391  grant. 

The  experimental  work  was  developed  in  the  Centre  for  Solar 
Energy  and  Hydrogen  Research  (ZSW),  Ulm,  Germany. 

The  authors  are  grateful  to  ZSW  people,  especially  to  Christoph 
Hartnig  to  coordinate  the  lab’s  work,  Juergen  Kaczerowski  for  the 
experimental  setup  and  Joachim  Scholta  for  his  suggestions  about 
the  measurement  system. 

Finally,  thanks  to  Michael  Danzer,  of  the  Ulm  University,  for  his 
valuable  discussions  and  feedback. 

References 

[1]  M.  Granovskii,  I.  Dincer,  M.A.  Rosen,  Int.  J.  Hydrogen  Energy  32  (8)  (2006) 
927-931. 

[2]  R.K.  Ahluwalia,  X.  Wang,  J.  Power  Sources  177  (2008)  167-176. 

[3]  D.M.  Bernardi,  M.W.  Verbrugge,  J.  Electrochem.  Soc.  139  (1992)  2477- 
2491. 

[4]  K.  Broka,  P.  Ekdunge,  J.  Appl.  Electrochem.  (1997)  27. 

[5]  V.  Gurau,  H.  Liu,  S.  Kakac,  AIChEJ.  44  (1998)  2410-2422. 

[6]  A.  Rowe,  X.  Li,J.  Power  Sources  102  (2001)  82-96. 

[7]  J.S.  Yi,  T.V.  Nguyen,  J.  Electrochem.  Soc.  145  (1998)  1149-1159. 

[8]  S.  Urn,  C.Y.  Wang,  J.  Power  Sources  125  (2004)  40-51. 

[9]  T.  Zhou,  H.  Liu,  J.  Power  Sources  138  (2004)  101-110. 

[10]  D,  Bevers,  M.  Wohr,  K.  Yasuda,  K.  Oguro,  J.  Appl,  Electrochem.  (1997)  27. 

[11]  M.  Ceraolo,  C.  Miulli,  A.  Pozio,  J.  Power  Sources  113  (2003)  131-144. 

[12]  B.R.  Silversen,  N.  Djalali,  J.  Power  Sources  141  (2005)  65-78. 

[13]  T.E.  Springer,  M.S.  Wilson,  J.  Electrochem.  Soc.  140  (1993)  3513-3526. 

[14]  T.E.  Springer,  Rainstrick,J.  Electrochem.  Soc.  136  (1989)  1594-1603. 

[15]  J.R.  Macdonald,  Impedance  Spectroscopy,  John  Wiley  8i  Sons,  Canada, 
1987. 

[16]  A.G.  Hombrados,  L.  Gonzalez,  M.A.  Rubio,  W.  Agila,  E.  Villanueva,  D. 
Guinea,  E.  Chimarro,  D.  Moreno,  J.R.  Jurado,  J.  Power  Sources  151  (2005) 
25-31. 

[17]  M.  Boillot,  C.  Bonnet,  N.  Jatroudakis,  P.  Carre,  S.  Didierjean,  F.  Lapicque,  Fuel 
Cells  1  (2006)31-37. 

[18]  N.  Wagner,  J.  Appl.  Electrochem.  32  (2002)  859-863. 

[19]  M.  Ciureanu,  H.  Wang,  J.  Electrochem.  Soc.  146  (11 )  (1999)  4031-4040. 

[20]  B.  Andreaus,  A.J.  McEvoy,  G.G.  Scherer,  Electrochem.  Acta  47  (2002) 
2223-2229. 

[21]  X.  Yuan.J.C.  Sun,  M.  Blanco,  H.  Wang.J.  Zhang,  D.P.  Wilkinson,  J.  Power  Sources 
161  (2006)  920-928. 

[22]  N.  Fouquet,  C.  Doulet,  C.  Nouillant,  G.  Dauphin-Tauguy,  B.  Ould-Bouamama,  J. 
Power  Sources  159  (2005)  905-913. 

[23]  A.  Carullo,  F.  Ferraris,  M.  Parvis,  A.  Vallan,  E.  Angelini,  P.  Spinelli,  IEEE  Trans. 
Instrum.  Meas.  49  (2)  (2000)  371-375. 

[24]  F.  Jurado,  J.  Power  Sources  129  (2004)  205-215. 

[25]  M.  Buchholz,  V.  Krebs,  Fuel  Cells  7  (5)  (2007)  392-401. 

[26]  YP. Yang,  F.CWang,  H.P.  Chang,  Y.-W.  Ma,  B.J.  Weng,  J.  Power  Sources 


1164(2007) 


MA.  Rubio  etal./ Journal  of  Power  Sources  183(2008 )  118-125 


125 


27]  M.A.  Rubio,  A.  Urquia,  S.  Dormido,  J.  Power  Sources  171  (2007)  670- 
677. 

28]  M.W.  Braun,  D.E.  Rivera,  A.  Stenman,  W.  Foslien,  C.  Hrenya,  Proceedings  of  the 
American  Control  Conference,  San  Diego,  1999. 

29]  L.  Ljung,  T.  Sbderstrom,  Theory  and  Practice  of  Recursive  Identification,  MIT 
Press,  1987. 

30]  Jer-Nan  Juang,  Applied  System  Identification,  New  Jersey,  1994. 

31  ]  J.  Larminie,  A.  Dicks,  Fuel  Cell  System  Explained,  John  Wiley  8;  Sons,  England, 
2000. 


[32]  R.F.  Mann, J.C.  Amphlett,  M.A.I.  Hooper,  H.M, Jensen,  B.A.  Peppley,  P.R.  Roberge, 
J.  Power  Sources  86  (2000)  173-180. 

[33]  M.  Usman  Iftikhar,  D.  Riu,  F.  Druart,  S.  Rosini,  Y.  Bultel,  N.  Retiere,  J.  Power 
Sources  160  (2006)  1170-1182. 

[34]  D.  Natarajan,  T.V.  Nguyen,  J.  Electrochem.  Soc.  148  (2001)  1324-1335. 

[35]  L.  Pisani,  G.  Murgia,  M.  Valentini,  B.  D’Aguanno,  J.  Power  Sources  108  (2002) 
192-203. 

[36]  I.  Sadli,  P.  Thounthong,  J.-P.  Martin,  S.  Rael,  B.  Davat,  J.  Power  Sources  156  (2006) 
119-125. 


